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OBJECTIVE 


The ultimate objective of this project is to develop a new technique which 
permits an accurate simulation of microstrip patch antennas or arrays with 
various feed, superstrate and/or substrate configurations residing in a re- 
cessed cavity whose aperture is planar, cylindrical or otherwise conformed to 
the substructure. The technique combines the finite element and boundary 
integral methods to formulate a system suitable for solution via the conju- 
gate gradient method in conjunction with the fast Fourier transform. The 
final code is intended to compute both scattering and radiation patterns of 
the structure with an affordable memory demand. 

PROGRESS SUMMARY 

This report contains four documents which represent our progress over the 
past six months. In summary, we continue to improve and develop our origi- 
nal code for other applications and we have begun the development of a new 
formulation and associated code for cylindrically conformal antennas /arrays. 
As stated in the included reports the array code developed last year was mod- 
ified to permit modeling of resistive cards, patch coatings and several new 
feed configurations, including striplines and aperture coupled feeds. Using 
these new capabilities we examined the RCS, input impedance, gain and res- 
onant frequency of several rectangular configurations using different loading 
and substrate/superstrate configurations. 

For the development of the new antenna code as applied to cylindrical 
arrays we require use of new discretization elements and expansion basis. 
One of the included reports in this binding describes the required elements 
and we are currently proceeding with the implementation of a formulation 
based on these elements. 

TRANSITIONS 

Over the next six months we shall concentrate on the development and val- 
idation of the code for antennas/arrays on a cylindrical surface. We shall 
also investigate efficient methodologies for extending this analysis to dou- 
bly curved conformal arrays. Further improvements and applications of the 
planar array code are also planned. 
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RCS Reduction of a Microstrip Patch 
Using Lumped Loads 

John L. Volakis and Angelos Alexanian 


Abstract 

In this report we consider the RCS of a single rectangular patch an- 
tenna in a recessed cavity. Using a previously developed finite element- 
boundary integral code, a study is performed on the patch’s RCS a & 
a function of frequency. To reduce the RCS of the patch at the reso- 
nant frequency, lumped (resistive) loads are placed at the edges of the 
patch. The effect of the lumped loads on the patch’s RCS and gain are 
examined and it is observed that the RCS and gain are reduced as the 
load value decreases, whereas the antenna’s bandwidth is increased. 
At resonance, the usual relation between RCS and gain is observed, 
but it is shown that this relation no longer holds at frequencies away 
from resonance. 
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Investigation’s Scope 

The purpose of this study is to (see figure 1): 

• Study the RCS vs. Gain of the microstrip antenna by placing lumped 
resistive loads at the side aperture of the patch 

• Determine the resonant frequency shift due to loading 

• Assess bandwidth enhancement due to loading 

This study was performed by using the University of Michigan program C AV- 
ITY3D (see Appendices). Other schemes for reducing the RCS of the patch 
and cavity edges using R-cards and dielectric coating are illustrated in figures 
1 and 2, but these have not yet been considered. 

Description of Investigation 

To assess the effect of lumped loads on the RCS and Gain of the microstrip 
patch we considered a 5cm x 3.4cm rectangular patch whose discretized form 
is shown in figure 3. This patch is loaded with four resistors placed roughly 
at the midpoint of each side as shown in figures 2 and 3. The thickness of 
the substrate/cavity was 0.17558cm and the relative dielectric constant was 
chosen to be e,. = 2.17. Also, for input impedance and gain computations 
the probe feed was placed at the bottom left corner of the patch as shown in 
figure 3. 

For the above patch, the resonant frequency was found to be /o = 
1. 94675GHz with the resistance of each load set at 15000 (essentially un- 
loaded). At this frequency we proceeded to evaluate the Gain and RCS of 
the patch as the resistance was changed from 15000 down to zero Ohms 
(shorted). The normalized values of these quantities as a function of load 
resistance are plotted in figure 4 and it is clear from this data that the RCS 
reduction is twice (in dB) that of the Gain (the RCS in figure 4 was computed 
with no load at the feed). This is to be expected since the RCS at resonance 
is given by [1] 

»(M) = -{(?(«,«}’ 1 

IT 

or 

10 log a(0, <j>) = 2 {lOlog[G'(0, <f>)\ + 101og[A/v'jr]} 
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where 10 log G' is the Gain in dB and includes any losses due to lumped or 
distributed loads. Let us assume that for a given load #1, the Gain is G\ 
and for another load #2, the Gain is G\. The corresponding cross sections 
evaluated at the same frequency are then given by 


10 log crj 


10 log a 2 


20 log G', + 10 log -p 
V* 


20 log G 2 + 10 log 




Subtracting these equations we obtain 

Aovib = 10(logcri — log cr 2 ) = 2[101ogGi — 101ogG 2 ] 


or 

Arras = 2[GidB — G 2 { ib] = 2AGdB 

which clearly shows that the dB RCS reduction is twice the dB gain reduc- 
tion. 

As expected, the input impedance becomes reactive with the addition 
of lumped loads and furthermore its real part is reduced in parallel with 
the lumped loads. Thus the patch’s resonant frequency changes with the 
value of the applied load. The specific values are illustrated in figure 5 
and the reason for the input resistance’s reduction is that the lumped loads 
placed between the patch edges and the ground are in a parallel arrangement 
with the impedance originally observed at the feed. However, we note that 
the radiation patterns are essentially unaffected by the placement of the 
loads (see figures 6(a) and 6(b)). More specifically, the radiation pattern 
corresponding to the 3000 loads is completely identical to that associated 
with the unloaded patch. We will thus choose to retune the patch with this 
load placed at its four side apertures (in actuality only the loads at two of 
the apertures are of significance since the field at the midpoint of the other 
two vanishes when the patch is at resonance). 

The resonant frequency (the frequency at which the input impedance is 
real) with the feed at the lower left corner (i = 1.25, y = 0.85 position) 
of the patch and with the 3000 loads at the patch edges was found to be 
/o = 1.97GHz. As given in figure 7 the input resistance at this frequency 
is approximately 103.50, and as shown in this figure, the input impedance 
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changes rather substantially as the feed position changes. Interestingly, when 
the feed is at (x = 2.5cm, y = 1.7cm), the input impedance is again real 
and equal to approximately 4911, a value which matches the characteristic 
impedance of common coaxial cables. From fig. 7, it is also observed that 
the gain is essentially unaffected by the feed’s position implying that the 
radiation pattern is also unaffected. 

With the feed at (x = 2.5cm, y = 1.7cm) and the 30011 loads at the patch 
edges, we next examined the patch VSWR and bandwidth. The results are 
given in figures 8 and 9, and the assumed source/line impedance was that 
at fo = 1.97GHz (i.e., Zo ~ 4911). It is seen that the patch bandwidth is 
increased substantially from fig. 8. Specifically, the bandwidth of the 30011 
loaded patch is approximately 0.13 (normalized frequency range within which 
VSWR< 3) whereas that of the unloaded patch is only 0.02 (more than six 
times smaller). The RCS of the 300H loaded patch in the vicinity of the 
resonant frequency is shown in figure 10 with or without a load of 4912 across 
the patch at the feed point. As expected, the RCS is maximum at the 
resonant frequency and then drops rather rapidly away from this frequency. 
When the 4912 load is placed at the feed to simulate the presence of the 
coax line, a 7dB reduction is observed at the resonant frequency. A more 
extensive frequency sweep of RCS for different loading configurations can be 
seen in figure 11. Here we are comparing an unloaded patch to one loaded 
with four 30012 resistors and to another loaded with four 30012 along with 
the 4912 at the feed (coax line simulation). Again we observe the expected 
RCS reduction at the points of resonance. 

Conclusions 

One can readily conclude from the above study that the RCS of a patch 
at resonance can be substantially reduced by placing lumped loads at its 
side edges. The reduction in RCS achieved in this manner is about twice 
(in dB) the corresponding reduction in gain. It was also noted that the 
resonant frequency is a function of the load’s value and thus the patch must 
be retuned for each load. Alternatively, it was shown that the feed location 
can be moved to again achieve resonance or simply to find an input resistance 
that matches that of the feeding line. 
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Future Study 

During the next month we shall investigate the RCS and gain of the patch 
when a dielectric coating is placed on its upper surface. Also the effect of 
a uniform R-card across the aperture will be studied. For the first case, a 
modification of the program is required. 
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Appendices 


Gain Calculations from Program Output Parameters 

By definition, 


Gain = G'(6, <f>) — lim 




p 


where E r is the radiated field as a function of the spherical angles 6 and <f>, 
Z 0 is the free space intrinsic impedance, r is the distance to the observation 
point and P m denotes the input power. For a probe feed, we have that 


Pin = \vr = ijy/ | 2 = \R^ 

since I = lAmp., in which R „ = Re { Z m } is the antenna input resistance. 
It then follows that 


G\9, <j>) 


7 D ^cm 


in which Acm denotes the wavelength in centimeters. The dB value of this 
expression is computed a 5 


GUM) = 101o glo 4r 


+ 101og 10 


' AL ' 

cm 

Z()Ri n 


The first term in this expression is that computed by the program CAV- 
ITY3D and is referred to in that program as the radiated power. The second 
term must therefore be added externally from the computed value of R m also 
given by the program. 


RCS Calculation from Program Output Parameters 

The RCS computed by the program CAVITY3D is given in dB above a square 
wavelength. To obtain the RCS in square meters the following formula must 
be applied externally: 


^dBSM — <7dBSW + 201og lo (A m ) 
where A m is the wavelength in meters. 
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RCS Reduction Studies 


Resistive cards 
Lumped loads 
Shape 

Feed mechanism 


Uniform Dielectric 


Substrate 


Tapered 
R Card 


Dielectric 
Coating on 
Patch 


Tapered 
r R Card 



Lumped 

Loads 



• Study RCS vs gain characteristics for various 
load treatments 

• Study resonant frequency shift due to loading 

• Study bandwidth enhancement as a function 
of loading 


Figure 1 


Other RCS reduction approaches 
for simple elements 



• Tapered R Card reduces RCS of terminations without 
effecting antenna parameters. 

• R Card strips act as lumped loads in reducing the RCS. 
They reduce RCS and Gain, but RCS is usually reduced 
much more. 

• Patch must be retuned after loading. 

• Radiation pattern is generally uneffected. 

• Surface wave effects are suppressed because of cavity 
enclosure. 


Goal: Given certain criteria with respect to Gain, bandwidth, 
and RCS at the design frequency, what should be the 
loading, feed mechanism and structural characteristics 
of the patch antenna? 


Figure 2 
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Figure 6(a) 
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Figure 6(b) 







Radiation Properties of Patch as Feed is Moved 
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Figure 7. The arrow on the patch above indicates the motion of the feed with 
initial and final positions at nodes (6,6>and (15,15) respectively. The frequency is 
1 .97GHz which happens to be the resonant frequency for four 300Q loads and the 
feed at node (6,6). With the feed at node (11,11) it is interesting to note that the 
input impedance becomes approximately 50Q (ie. 48.92-jl.13 Q). 


VSWR 



Frequency in GHz 


Figure 8. VSWR of loaded and unloaded patch vs frequency 
Computation of VSWR : 

VSWR = ( 1 +iri)/( 1 -ID) where V = (Zin-Zinr)/(Zin+Zinr) 

Zinr= 48.92-jl.13 £1 is the input impedance at the resonant frequency (1.97GHz) 
for four 300Q loads and feed at node position (11,1 1). Note that at this frequency 
VSWR=1. 

Zin = input impedance for the same configuration as above, with the exception that 
the frequency varies. 





Frequency in GHz 

Figure 10. RCS as a function of frequency 
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Figure 11. RCS as a function of frequency 




Typical input file 


7.5 5.1 .17558 
31 31 2 
1 

(2.17,0) (l.,0.) 
1.97 1.97 1 
2 

0 0 1 
0 0 1 
.01 800 
2 
1 

6 6 
20 20 
1 
2 
0 
1 

11 11 
1 0 
1 
4 

15 6 
(300., 0.) 

1 

25 15 
(300. ,0.) 

1 

15 25 
(300. ,0.) 

1 

6 15 
(300., 0.) 

1 

2 


cavity size in cm (x,y,z) 
number of nodes (x,y,z) 

- uniform dielectric 
substrate dielectric - er pr 
frequency min, max, increment 

- radiated intensity calculation 
theta min, max, increment 
phi min, max, increment 
tolerance, number of iterations 

- monitor cnvergence l=yes 2=no 
number of patch 

lower left comer of patch node(x,y) 
number of elements on patch (x,y) 
layer in which patch resides 

- patch display l=yes 2=no 
number of short circuit pins 
number of feeds 

feed position node(x,y) 
amplitude and phase of current source 
layer in which feed is embedded 
number of loads 
node(x,y) load 1 position 
value of load 1 (complex) in ft 
layer in which load 1 is embedded 
node(x,y) load 2 position 
value of load 2 (complex) in ft 
layer in which load 2 is embedded 
node(x,y) load 3 position 
value of load 3 (complex) in ft 
layer in which load 3 is embedded 
node(x,y) load 4 position 
value of load 4 (complex) in ft 
layer in which load 4 is embedded 

- resistive card l=yes 2=no 


Figure 12 




Actual element dimensions for the following patch and cavity 
configurations. 
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Patch dimensions : x-5.0cm. y-3.4cm 


Figure 13. This is the exact patch-cavity configuration resulting from the input 
file given in figure 12. 
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Figure 14. Extra grid for future designs. 
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A scheme to alter the resonant frequency 
of the microstrip patch antenna 

J. L. Volakis and J. M. Jin 
December 20, 1991 



Abstract 

Simple schemes are presented for altering the resonant frequency 
of the rectangular patch antenna without a need to change its size. 

In particular, by placing a perturbance below the patch it is shown 
that as much as 20 percent increase and 30 percent decrease from 
the resonant frequency of the unperturbed patch can be achieved. 

The specific configurations considered in this letter include a cavity- 
backed, aperture-backed and a protrusion-backed patch, and for each 
case design curves are presented. 

Introduction 

In many cases one is required to design a certain antenna or array which 
does not exceed given physical dimensions. Since the size of the antenna is 
typically determined by its operational frequency, restrictions on the physical 
dimensions of the antenna can present a challenge to the antenna designer. 
To facilitate the design process, in this letter we consider simple schemes to 
rather substantially lower or increase the operational frequency of the rect- 
angular microstrip patch antenna without altering its aperture size. These 
include the placement of some perturbance (protrusion, depression or slot) in 
the cavity region below patch, a scheme well known to designers for control- 
ling the cavity’s resonant frequency [1], [2] or other antenna parameters [3]. 
However, because of a lack of available CAD tools to treat such pertur- 
bances, the approach does not appear to have been explored for altering the 
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frequency of the highly resonant microstrip antennas. Moreover, traditional 
applications of these frequency control methods dealt with small shifts in the 
resonant frequency. 

Recently, a new analytical technique was developed for a characteriza- 
tion of microstrip antennas or arrays residing in a cavity that is recessed in 
a ground plane [4]. The technique is based on a hybrid version of the finite 
element method (FEM) where the mesh is terminated at the aperture of 
the patch via the exact boundary integral equation. Because of the FEM’s 
geometrical adaptability perturbances or irregularities in the cavity’s geom- 
etry, and material inhomogeneities in the substrate can be modeled without 
difficulty. Below we specifically consider the effect of three different cav- 
ity configurations on the resonant behavior of a given patch residing at the 
aperture of the cavity. 

Antenna Geometry 

Consider the reference geometry illustrated in Fig. 1. The configuration 
consists of a square 4.45cm x 4.45cm patch residing in a square cavity whose 
aperture size is 6.625cm x 6.625cm and depth is .05cm. The cavity is filled 
with a dielectric having a relative permittivity of e r = 2.17 and we remark 
that its finite size has negligible effect on the resonance and input impedance 
of the patch. 

With the goal of altering the resonant behavior of this patch we consider 
three modifications of the cavity. The xz and yz cross sections of these are 
illustrated in Fig. 2(a)-(c). They include a rectangular depression at the 
center of the original cavity, a second cavity which is aperture coupled to the 
original and a third configuration having a rectangular protrusion below the 
patch. In all cases the patch and aperture size is unchanged from the original 
antenna shown in Fig. 1. 

The analysis of all four patch antennas depicted in Figs. 1 and 2 was ac- 
complished via the hybrid finite element-boundary integral (FE-BI) method 
described in [4]. In accordance with this method the cavity region below the 
patch is formulated via the finite element method, which is readily adaptable 
to any cavity configuration with an appropriate choice of element shapes. For 
the particular cavity geometries at hand we employed rectangular bricks with 
edge-based shape functions which have zero divergence and are appropriate 
for modeling metallic corners. Because the FEM solves the wave equation, it 
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yields a solution subject to the boundary conditions at the cavity walls and 
its aperture. For our case, the tangential field is set to zero at the interior 
cavity walls and at the aperture we enforce the field continuity between the 
interior and exterior fields. The last fields are expressed in terms of the ra- 
diation integral with the sources set equal to the tangential aperture fields. 
Consequently, the continuity condition on the aperture fields results in an 
integral equation which can be alternatively thought as a global mesh termi- 
nation condition. In contrast to the usual absorbing boundary condition(s) 
which must also be enforced far from the radiating element, this approach 
is exact but involves a partially full submatrix. However, by solving the 
resulting (sparse) system via the biconjugate gradient method the need to 
generate the matrix is eliminated thus retaining the O(N) storage require- 
ment, which is typical in FEM solutions. As noted in (4) we have found this 
solution approach to be highly efficient with respect to the system size and 
convergence rate. 


Numerical Results 

Plots of the resonant frequency associated with the configurations in Figs. 
2(a) and 2(b) to be referred to as cavity-backed or aperture-backed antennas, 
respectively, are given in Figs. 3 and 4. In Fig. 3, we depict the resonant 
frequency as a function of the second (bottom) cavity depth with the aperture 
at the interface of the two cavities fixed at 3.56cm x 3.56cm. The plot in 
Fig. 4 illustrate the effect of the aperture size on the resonant frequency with 
the bottom cavity depth kept at 0.3cm. From both figures, it is seen that 
by changing the aperture or depth of the bottom cavity one can decrease 
the patch’s resonant frequency by as much as 30 percent. Consequently, 
one may construct a 30 percent smaller patch from the conventional one 
with the introduction of the second cavity. We have also verified that the 
RCS of the new antenna structure is nearly the same as that of the original, 
and thus the new design results in an overall reduction of the antenna’s 
radar cross section (RCS). It should be noted though that when the aperture 
of the second cavity becomes larger than the patch size, there is minimal 
frequency shift and this is clearly shown in Fig. 4. The shift in the resonant 
frequency can be understood by resorting to a transmission line model of 
the patch antenna. In this case the inserted aperture below the patch can 
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be represented by an equivalent reactance placed between the admittances 
representing the patch terminations. 

The resonant frequency of the third antenna (protrusion-backed) depicted 
in Fig. 2(c) is plotted in Fig. 5 as a function of the protrusion’s height. For 
these calculations the xy cross section of the protrusion was fixed to be 3.56cm 
x 3.56cm. Not surprisingly, the protrusion increases the resonant frequency 
of the original patch almost linearly and up to 20 percent as the protrusion 
height is increased. As expected, when the protrusion’s cross section extends 
beyond the patch, there is minimal shift in the resonant frequency of the 
reference patch having a substrate thickness equal to the distance between 
the patch and protrusion. 

Unfortunately, the bandwidth of the cavity-backed and aperture-backed 
configurations was not altered although some differences were observed in 
the actual values of the input impedance as shown in Fig. 6. To increase 
the patch’s bandwidth one could, of course, resort to standard approaches 
such as those reported in [5], [6]. Some bandwidth increase is observed for 
the protrusion- backed antenna but this is attributed to the larger substrate 
thickness below the patch’s edges. 
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Figure 1 . Reference geometry of the microatrip patch in a metal-backed rectangular 
cavity. 





Figure 2. Modified cavity configuration* r- (a) cavity- backed, (b) aperture- backed, 
(c) protrusion-backed. 


x> *> 



2300 . 
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Figure 3. Resonant frequency of the aperture- and cavity-backed patch configura- 
tions as a function of the bottom cavity depth d . The original cavity depth 
is retained at 0.05cm and the entire cavity region is filled with a dielectric 
having e r = 2.17. The aperture of the second cavity is 3.56cm x 3.56cm. 
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Figure 4. Resonant frequency of the aperture-backed patch as a function of the 
square aperture dimension w. The depth of the bottom cavity was kept 
at d = 0.3cm and all other parameters are the same with those used for 
Figure 3. 





For (c), the feed was placed at the diagonal of the patch and a distance 
of 1.259cm from the nearest corner. 
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Figure 5. Resonant frequency of the protrusion-backed patch as a function of the 
protrusion’s height h. The protrusion’s xy cross section was kept at 3.56cm 
x 3.56cm. 
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A Collection of Edge-Based Elements 

Leo C. Kempel and John L. Volakis 

Radiation Laboratory . \ ) 

University of Michigan / 

January 5, 1992 


Abstract 

Edge-based elements have proved useful in solving electromag- 
netic problems since they are divergenceless. Previous authors have 
presented several two-dimensional and three-dimensional elements. 
Herein, we present four types of elements which are suitable for mod- 
eling several types of three-dimensional geometries. Distorted brick 
and triangular prism elements are given in cartesian coordinates as 
well as the specialized cylindrical shell and pie-shaped prism elements 
which are suitable for problems best described in polar cylindrical 
coordinates. 



Contents 

1 Introduction 

2 Distorted Brick Element 

3 Distorted Triangular Prism 

4 Cylindrical Shell Edge-Based Elements 

5 Pie-Shaped Prism 
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Conclusions 


where E is the electric field and G is a Green’s function. For the remainder 
of this report, we suppress the dependence of F on the unprimed coordinates 
since there will be no operations carried out in these coordinates. Further- 
more, we express (1) in terms of its cartesian components 

7 = xl x + yl y + zl x (3) 


where 

h = C • J v E(x',y',z')dx' dy' dz' (4) 

and C € { x,y,z }. 

To evaluate (4), we apply the change of variables theorem to map the 
original brick domain into the unit cube as shown in figure lb. 


I( = C • / / / F(x'(u,v,w),y'(u,v,w),z'(u,v,w))\J(u,v,w)\dudvdw (5) 
J 0 J 0 J 0 


which is accomplished by allowing 
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x'(u, V, w) 

= '52N i {u,v,w)x , i 

( 6 ) 


t=l 
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y'(UyVyW) 

= Ni(u,v,w) y[ 

1=1 

( 7 ) 

Z'(UyVyW) 

= ^2 Ni(u,v,w) z\ 
1 = 1 

( 8 ) 


with (x^y'^z'i) denoting the ith node coordinates as shown in figure la. The 
shape functions are given by 


Ni(u,v,w) 
N 2 (u,v,w) 
Nj'u,v, w) 
N 4 {U, Vy W ) 


(1 — u)(l — v)(l — w) 
u(l — n)(l — w ) 
uv(l — w) 

(1 — u)n(l — w) 


N s (u,v,w) 
N 6 (u,v,w) 
Nj(u, v, w ) 
Ns{u,v,w) 


and the Jacobian is of the form 


(1 — u)(l — v)w 
u(l — v)w 
uv w 

(1 — u)vw 

( 9 ) 


_ MtyW _ di i'_dx!_dE_ 

du dv dw du dw dv du dv dw 
dy' dx' dz' dz 1 dx' dy' dz 1 dx 1 dy ' 

^ du dw dv du dv dw du dw dv 


( 10 ) 
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1 Introduction 

Edge-based expansion functions have become popular in finite element and 
moment method techniques applied to the solution of electromagnetic scat- 
tering problems. This is due to the fact that the expansion functions are 
divergenceless within the element domain while any fictitious charges formed 
at the edges of the element will be canceled by another edge resulting in a 
volume which is divergenceless. In contrast, nodal expansion functions are 
not divergenceless and are thus less accurate. 

Tanner and Peterson [1] presented triangular and tetrahedral divergence- 
less vector elements while Jin and Volakis [2] applied the corresponding rect- 
angular brick elements to slot scattering problems. Kameari [3] presented 
distorted linear and quadratic brick elements in the context of transient 
quasi-magnetostatic problems. 

In this report we derive edge-based expansion functions for a variety of 
three-dimensional elements. These expansion functions are characterized by- 
unit value along one edge of the element while vanishing along all the re- 
maining edges which are parallel to the unit edge. Distorted bricks, distorted 
triangular prisms, cylindrical shells and pie-shaped prisms will be presented. 
The first two of these are suited for completely general inhomogeneous do- 
mains whereas the last two are more appropriate for domains terminated at 
cylindrical boundaries. 

2 Distorted Brick Element 

The first element considered is the distorted brick which is shown in figure la. 
The strategy used herein is to first employ the change of variable theorem to 
convert the integration domain from the original distorted brick into a cube. 
Once this is accomplished, the vector-valued function may be expanded in 
terms of edge-based expansion functions as was done by Jin and Volakis[2]. 
Typically, one wishes to integrate some vector-valued function 

I(x,y,z ) = J v F(x',y',z',x,y,z)dx' dy' dz' (1) 

over the domain of the element where V is the volume of the brick and the 
vector valued integrand is typically of the form 

F(x\ y\ z') = £(*', y', z') G(x\ y\ z\ x,y,z) (2) 
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where E is the electric field and G is a Green’s function. For the remainder 
of this report, we suppress the dependence of F on the unprimed coordinates 
since there will be no operations carried out in these coordinates. Further- 
more, we express (1) in terms of its cartesian components 

I = xl x + yl y + zl z ( 3 ) 

where 

■ J v F(x',y',z')dx' dy' dz' ( 4 ) 

and C € {x,y,z}. 

To evaluate (4), we apply the change of variables theorem to map the 
original brick domain into the unit cube as shown in figure lb. 

h-C‘f f f F{x'{u,v,w),y'(u,v,w),z'(u,v,w))\J(u,v,w)\dudvdw (5) 

JO JO JO 

which is accomplished by allowing 

s 

x'(u,v,w) = ^ Ni(u } v, w) x\ (6) 

t=i 

8 

y'(u,v,w) = ^ Ni(u, v, w) y[ (7) 

1=1 
8 

z'(u,v,w) = ^2 N,(u, v, w) z\ (8) 

t=l 

with (x^, y,, z x ) denoting the ith node coordinates as shown in figure la. The 
shape functions cire given by 


Ni(u,v,w) = (1 — u)(l — u)(l — w) Ns(u,v x w) 

N 2 (u,v,w) = u(l - t>)(l - w) jV 6 (u, v, w) 

N 3 (u,v,w) = uu(l-tu) N 7 (u,v,w ) 

N 4 (u,v,w) = (1 — t*)v(l - w) Ns(u,v,w) 

and the Jacobian is of the form 


(1 — u)(l — v)w 
u(l — v)w 
uvw 

(1 — u)vw 

( 9 ) 


dx' dy' dz' dx' dy' dz' dy' dx' dz' 

du dv dw du dw dv du dv dw 

dy' dx' dz' dz' dx' dy' dz' dx' dy' 

du dw dv du dv dw du dw dv 
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The partial derivatives are computed from 

c?C dNi . . ( 1 1 ) 

with C € { *, y, * } and i/ € { «, t>, to }• From (9) the partial derivatives of the 
shape functions are found to be 


du 

iV 2 (u. v, w) 


du 

diy 3 (u, v, w ) 

du 

dN 4 {u,v.w) __ 

du 

dNiju , v, w ) _ 

dv 

dNi{u,v,w) _ 

dv 

dNz{u,v, w) _ 

dv 

dN^u, v, w) _ 

dv 

dNi{u,v,w) __ 
d w 

dN^ju, v, w) _ 

dw 

dN^ju^w^ __ 

dw 

dN 4 {u,v,w) 

dw 


-(1 - u)(l - to) 

dN s {u,v,w) 

du 

= - (1 —v)w 

(1 — v) (1 w) 

dN e (u, v, w) 

du 

= (1 — v) w 

v (1 — w) 

dN 7 {u,v,w) 

du 

= VW 

— v (1 — w) 

dN&{u,v,w) 

du 

— — VW 

-(1 - u)(l - w) 

dN s {u,v,w) 

dv 

— — (1 — ti) w 

-u (1 — w) 

dN 6 (u,v,w) 

dv 

= —inv 

u (1 — w) 

dN 7 (u,v,w) 

dv 

— uw 

(1 - u) (1 - w) 

dN&{u,v,w) 

dv 

3 

"S' 

1 

rH 

II 

-(1 -u)(l -v) 

dN S {u,V,w) 

dw 

= (l-u)(l-v) 

—u (1 — u) 

dN e (u, v , w) 

dw 

= u (1 — t>) 

—uv 

dN 7 (u,v,w) 

dw 

= uv 

- (1 - u) v 

dNa{u,v,w) 

dw 

= (1 — u) V 
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Now that we have transformed the integration domain from the original 
distorted brick into the unit cube, we may expand the vector-valued integral 
in terms of an edge- based expansion of the form 


12 

F(x'(u,v,u>)V(u,^t<Oi *'(“»”»«>)) = 5ZiV*(u,u,u;)0j (12) 

j= i 


where 


A 7 K 

v, 

, id) = 

(1- 

”)(1 A - 

w)t 12 

NJ(u, 

u, 

w) = 

u(l - 

W)t 58 

NZ(u, 


id) = 

d(1 - 

- w )<43 


A r |(u, 

v, 

w) = 

uwt 67 



U, 

id) = 

(1- 

V)wt 56 


Ni(u, 

v, 

id) = 

(1 — U 

)(i -v)i 

A r 4 e (u, 

V, 

w) — 

VU't S 7 


NJ 0 (v, 

v, 

id) = 

u(l — 

v)t 26 

JV 5 e (u, 

V, 

id) = 

(1- 

«)(!“ 

w)t 14 


V, 

id) = 

(1 — U 

>*48 

Ni(u, 

V, 

id) = 

(1- 

u)wi 2 3 


jVj 2 (u, 

v, 

id) = 

UVt 37 



(13) 

in which are the unknown expansion coefficients and t XJ denote the ele- 
ment’s unit edge vectors given by 





(14) 


with r } being the position vectors in the original domain. 
We may now express the cartesian components of (3) as 


A 1 ^ y 1 yl yl _ 

= C • / / / jV e (u, u, tu)| J(u, u, ty)| <fu dudu; (15) 

;_1 Jo Jo Jo J 

recognizing that ifrj represents the average component of the vector-valued 
function associated with the jth edge. 

It should be noted that the brick cannot be too distorted. This is due to 
the change of variables theorem. In order to assure a one-to-one mapping for 
(x, y, z ) to (u, v, w ), the sign of the Jacobian may not change. This condition 
will preclude elements that are twisted or extremely narrow and diamond 
shaped. The typical element should therefore not deviate much from the 
rectangular brick. 
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3 Distorted Triangular Prism 

A triangular prism element such as the one shown in figure ‘2a may be de- 
veloped in exactly the same manner as the brick. In this case, we employ 
a mapping from the original volume V to the volume shown in figure 2b 
followed by an edge-based expansion of the vector-valued integrand. The 
cartesian component integral may be expressed as 


/=(./ [ [ F(x'(u,v,w),y'{u,v,w),z\u,v,w))\J{u,v,w)\dvdudw 

Jo Jo Jo 

where the coordinate mapping is given by 

6 

x'(u,v,w) = £ M(u, v, to) xj 

1=1 
6 

y'(u,v,w) = ^,Ni(u,v,w)y'i 

1 = 1 
6 

z'(u,v,w) = ^2Ni(u,v,w)z- 


(17) 

(18) 
(19) 


t = l 


with the node coordinates (x',yj,z') shown in figure 2a and the shape func- 
tions are given by 


Ni{u,v,w) = (1 - u — v)(l — to) N 4 (u,v,w) 
N 2 (u,v,w) = u(l — w) N s (u,v,w) 

N 3 {u,v,w) = v(l-ttf) N 6 (u,v,w) 


The Jacobian is again given by (10) with 


d£ * m 

dv rg \ dv 


( u,v,w)(i 


(1 — u — v)w 

uw 

vw 

( 20 ) 


(21) 
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and the non-zero shape partial derivatives of iV t are 


<9A r i 

n n dN ' 

= (1 w )= a 

c iv 

dN, 

u — v) — - 

dN 4 

du 


dw 

dN 2 

II 

» — 1 

1 

II 

dN? 

dN s 


du 

dw u ~ 

dw 


dN 4 

II 

1 

6 

II 

1 <*> 

II 

I 

II 

£ 

CO 

dN e 


du 

dN s 

dv 

dN 6 

aw 

dw 



^iV 5 (7iVg 


The vector valued function in (16) is expanded in terms of edge-based func- 
tions 

9 

F{x'{u, v , u>), y\u, v, w), z'(u, v, w)) = £ iV?(u, v, w ) (23) 

i=i 

where the vector expansion functions are 


N‘(u, v, w ) 
N%(u, v, tt;) 
N£(ti,v,w) 
N‘(u,v,w) 
N*(u,v,w) 
N£(u,v,w) 
N‘(u,v,w) 
iV|(u, v,w) 
N£(u,v,w) 


[(£l + ^2)^12 + ^2^13](1 — U>) 

— V^[^3^12 — ^2^13](1 — **>) 

— [&}<12 + (£l + ^3)<13](1 — w) 
[(6 + ^2)^45 + £2*4e]u> 

— \/2[^3^4S + ^2^4e]^ 

— [^3^45 + (£1 + £3)£w]u> 

u(l - v)t 2 & 

(1 - u)(l - v)tu 
(1 - 12)^36 


(24) 


where the area coordinates in (24) are given for the new integration domain 


— 2ysjj(sji <2 n )(s n f*n)(^n Cn) 

3n = l(a„ + 6 n + Cn) (25) 

« € {1,2,3} 
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and 


ai = y j(u ~ l) 2 + 

bi = + (u - l ) 2 

ci = n/2 

We are now able to exp 


a 2 = ^/u 2 + (v - l) 2 

i 2 = \/u 2 + u 2 

c 2 = 1 

(16) cis 


a 3 = \/u 2 + u 2 

63 = ^/(u - l) 2 + t- 2 

C 3 = 1 

(26) 


J N‘{u,v,w)\J(u,v,w)\dvdudw (27) 

J = 1 0 0 

realizing that the same limited distortion constraint exists for the triangular 
prism as was the case for the brick. 


4 Cylindrical Shell Edge-Based Elements 

Although the brick and triangular prism elements given previously are quite 
general, we find it advantageous to develop an element which is designed 
for cylindrical problems. The resulting formulation will be quite compact 
compared to the brick for example since a specific geometry is assumed. This 
element is the polar cylindrical coordinate equivalent of the brick introduced 
by Jin and Volakis [2] and is shown in figure 3. 

For cylindrical problems, we wish to evaluate the integral 

7= f f f F{p, 4>, z)p dp d<f>' dz' (28) 

J Z\ 1 J P\ 

where the vector- valued function F(p\ ft, z f ) is analogous to (2). We do not 
employ any change of variable for this element so we may directly express 
the vector- valued function in terms of an edge-based expansion 
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The shape functions are given by 


' p\ - p' 2 


V - z s 

.pI-pI. 


Zl - z 5 


<t> 


(29) 
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where (p x ,0 x 
is given by 


A r |(p\0V) = 

Nl{p\0\z') = 
N' 4 (p',0\z') = 

n;( p \0\z') = 

iV|(p\0',*') = 

N' 7 {p\4>',z') = 
W,0',z') = 

Nl{p',0\z') = 

n; 0 (p\ 0 ', z') = 

$ e n{p\4>',z') = 

N' 12 (p\4>\z') = 


1 - 


102 - <j/| 


102 ~ 01 1 


‘5 


*1 ~ z§ 


V’ ~ P? 


P\ ~ Pi 


z “ z 5 


1 - 


Z 1 ~ z 5 

10' - 01 1 


<P 


1 02 — <^1 1 J 



»D 

1 

l 


Z \ “ -^5 J 


P\ ~ P' 


2 — n 2 
4 Pi 

1 02 ~ 0'\ 



z' - z l 


, z l - 25. 


1 - 


1 02 ~ 01 1 J 


P' 2 ~ P\ 


4 ~ Pi 

f ~0l\ 



1 

1 

N 


.Zs-Zi\ 


0 

Z_ - Zl 
L 2 5 ~Zl\ 

0 


1 - 



z' - Z X 


z 5 — Z\ 


P4~ P 


Lp4 - Pi 
P 4 ~ P' 


LP4 — PlJ 

P' ~ Pi 


102 ~ 0'| 
102 “ 01 1 J 
10' -01 1 


P 4 ~ Pi 
P' - Pi 


1^2 — 1 J 

10' -01 1 


P 4 ~ Pi 


[102 — 01 1 
102 ~0'| 


102 — 01 1 J 


(30) 


t ,Zi) denote the node coordinates as shown in figure 3 and (28) 
r*b r<h fPi _ 


' = e rcr *'• 2< ) />' 


(31) 
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5 Pie-Shaped Prism 

The final element considered is the complement to the cylindrical shell ele- 
ment, the pie-shaped element which is shown in figure 4. The development 
of this element closely parallels the previous element. 

The vector- valued integral is given by 


/= rrrW'Msvw 

J Z\ J <t>2 J P\ 

The vector-valued function is expanded in terms of edge vectors 

F{p',4>',z') = ^N‘{p\4>\z')fa 
i = l 


and the edge-based expansion functions are 


N{{p,fa, z ') 

Nl(p\4>',z') 

Nl{p',<t>',z') 

Nl(p',4>\z') 

Nl{p',4>',z') 

N t 7 {p\4>\ z ') 

Nl{p',<t>',z) 

NZ(p',<t>',z) 


T - 



1 

\fa ~ 

*1 

fall 

\ Z *~ 

l*. 

-'1 

4> 


l Z * ~ Z l. 


\fa - 4>'\ 


\fa — fa\. 

| 1 j i z 4- 



(32) 


(33) 


(34) 
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where (/?,,©,, r,) are the nodal coordinates shown in figure 4. With this 
expansion. (32) may be written 

/=!>, rrr N’(i>‘A',z')p'dp'd4,'dz' ^ 

j— i J *\ J Pi 

6 Conclusions 

Four different linear volume elements were presented: distorted brick, dis- 
torted triangular prism, cylindrical shell, and pie-shaped prism. All of these 
elements are used to evaluate a vector- valued integral over the volume of the 
element by employing edge-based vector expansion functions. These func- 
tions were explicitly given for all elements. In addition, the shape functions 
and Jacobian associated with the mapping of the arbitrary brick or triangular 
prism to a known prism are presented. Since these elements are divergence- 
less, we expect them to useful for simulating various geometries via the finite 
element method or the Method of Moments. 
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Abstract 

A formulation is presented for modeling a resistive card in the 
context of the finite element method. The appropriate variational 
function is derived and for validation purposes results are presented 
for the scattering by a metal-backed cavity loaded with a resistive 
card. 
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1 Introduction 


A resistive card is an infinitesimally thin sheet of material which allows par- 
tial penetration of the electromagnetic field. Thin dielectric layers and very 
thin conductors whose thickness is less than the skin depth are examples of 
materials which can be modelled by resistive cards or sheets. Resistive cards 
are often used for radar cross section and RF power penetration control and 
as a result they have been studied extensively. Such studies have generally 
been done in the context of high frequency [1,2] and integral equation solu- 
tions [3,4], but to date the treatment of resistive cards within the context of 
the finite element method (FEM) has not been considered. Over the past 
few years, FEM has been applied to a variety of electromagnetic applications 
and it is thus important to incorporate the modeling of resistive cards in the 
FEM. In this paper we propose an FEM formulation which accounts for the 
presence of resistive sheets. To validate this formulation, results based on a 
physical modeling of the resistive sheet are also presented. In this case, the 
resistive sheet is equivalently replaced by a thin dielectric layer. The mod- 
eling of such a layer in the usual manner leads to larger and consequently 
inefficient linear systems, which is the primary reason for resorting to a math- 
ematical modeling of the resistive sheet. Results based on the mathematical 
and physical modeling are presented in connection with the scattering by a 
metal-backed cavity in a ground plane and these are used to validate the 
proposed mathematical model. 


2 Formulation 

A resistive sheet is characterized by its resistivity R which is measured in 
Ohms per square. Mathematically, it satisfies the boundary condition [5 

n x (h x E) = —Rn x (H + - H - ) (1) 

where denotes the magnetic field above and below the sheet, E is the 
electric field and its tangential component is continuous across the sheet, and 
n denotes the unit normal to the sheet pointing in the upward direction ( + 
side). To a first order, this boundary condition can be used to simulate the 
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presence of a thin dielectric layer by setting [4,5] 


R = 


Zo 


jk 0 (c r - 1) t 


In this, t is the thickness of the layer, Z Q and k 0 denote the free space intrin- 
sic impedance and wave number, respectively, and t r is the relative dielectric 
constant of the layer. Alternatively, a resistive sheet may be equivalently re- 
placed by a thin dielectric layer having thickness t and a relative permittivity 
of 


€r = 1 ~ 


jZ 0 

k 0 tR 


( 2 ) 


Generally, the accuracy of this simulation increases as the thickness t is 
decreased. Typically, t should not exceed one-tenth of the wavelength in the 
material. 

Let us now consider a finite element solution of the fields within a volume 
V subject to a given excitation. The volume consists of some inhomogeneous 
dielectric having relative permittivity and permeability t r and /x r , respec- 
tively, and we shall also assume that resistive cards may be embedded within 
the dielectric (see Figure 1). In accordance with the finite element method, 
the volume is subdivided into M smaller volume elements and in this case 
we require that the resistive cards are tangential to the boundary surface 
of these elements. A weak solution of the fields within the volume can be 
obtained by extremizing the functional 


M 


e=l 


( 3 ) 


F e 


///v.t (VxE) ' (V><E| -*- £ ' E ' E 

4 -jk o Z 0 (j^f E • (H x h c ) ds 


dv 


( 4 ) 


with respect to the electric field E including that implied in H. In this 
expression, V e is the volume of the element which is enclosed by the surface 
S e and h c denotes the outward normal to S € . 
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Generally, for a dielectric volume not enclosing resistive sheets or other 
current sheets, the contributions of the surface integrals in (3)-(4) vanish 
everywhere except when S e coincides with the outer surface S a of the volume 
V. This is a consequence of field continuity across the elements, but if a 
portion of the element’s surface coincides with a resistive sheet, then the 
surface integral in (4) does not vanish since the magnetic field is discontinuous 
as described in (1). Let us for example consider the surface S re which borders 
the eth and (e + l)th elements, and is coincident with a resistive sheet of 
resistivity R. Then the contribution from this surface to the surface integral 

in (4) is 

jk o Z 0 II E • (H- x n re ) ds 
J J Sre 

from the eth element and 


—jk 0 Z 0 J J s E • (H + x n re ) ds 


from the (e + l)th element with n Te pointing from the eth to the (e + l)th 
element. Combining these two integrals and employing (1), it follows that 
the contribution of the surface S rc to the surface integral in (4) is 


j ko 


./ 1 JE • (H* - H-) x ft,,] 


ds 


= jk 0 Z 0 J J $ J (« re x E) • (n re x E) ds (5) 

Consequently, the functional F may be rewritten as 

f - s///,t (V x E) • (V x E) - k 2 0 e r E • e| dv 
+jk 0 Z o J J s j (n, x E) • (h T x E) ds 

+jk 0 Z 0 (fj> s E-(Hxh 0 )ds ( 6 ) 

in which S T denotes the surface occupied by the resistive sheet and S 0 is the 
outer surface enclosing the volume V. As usual, n r is the outward normal to 
S r and h 0 is correspondingly the outward normal to S a . If S T borders the 
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outer surface of the volume V, then S 0 should be considered to be just over 
the exterior side of S r (i.e., S 0 always encloses S r )- 

Having derived the explicit form of the functional F, we may now expand 
the element field using the standard linear shape functions. If the sources 
are within V then F should be modified to read 


dv 


J J / v ,{^-( VxE ) ( VxE ) - k 2 a e r E • E j dv 

+ J J J v E • jk o Z 0 J mt - V x ^M* nt j 
+jtc 0 Z 0 J ^(n r x E) • (n T x E )ds 

+jk 0 Z 0 (j£s E • (H x h 0 )ds 


( 7 ) 


where (J m, ,M' n ‘) denote the impressed sources internal to V. Then, upon 
setting the first order variation of F to zero, we can obtain a system of 
equations for the solution of the interior and boundary electric fields. For a 
unique solution of this system we must, however, specify a relation between 
the tangential electric and magnetic field which appear in the surface integral 
over S 0 - If we assume that the subject volume is that occupied by the metal- 
backed cavity recessed in a ground plane, as shown in figure 2, then S Q reduces 
to the aperture area of that cavity. By invoking image theory, the magnetic 
field on the aperture can then be expressed as 


H = H ,nc + H 


ref 


- 2An // s .[( !+ 4 W ) G “ (r ' r,) 


[E(r) x z]ds (8) 


where S a denotes the aperture surface, G 0 (r,r') is the free space scalar 
Green’s function, r specifies the observation point located on S a and I = 
xi + yy + zz is the unit dyad. Also. H m<: denotes the magnetic field gen- 
erated by sources in the free space and H re ^ is the corresponding reflected 
field when the cavity’s aperture is shorted. Substituting (8) into (7) gives 
a functional only in terms of the electric field. The system obtained from 
this functional will be partly sparse and partly full. In particular the volume 
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integrals and that over S r in (7) lead to a sparse submatrix involving the 
interior Helds of the cavity. However, in view of (8) the last surface integral 
of (7) over S a (or S a ) renders a full Toeplitz submatrix since the boundary 
integral is convolutional. Consequently, by resorting to an iterative solu- 
tion such as the conjugate or biconjugate gradient method in conjunction 
with the fast Fourier transform, the need to generate the Toeplitz matrix is 
eliminated thus maintaining the 0(n) storage requirement, characteristic of 
finite element solutions. The details pertaining to this implementation are 
discussed in [6] - [8]. In the next section we only present some results aimed 
at evaluating the accuracy of the proposed resistive sheet model. 


3 Numerical Results 

Let us consider the metal-backed rectangular cavity illustrated in Figure 2. 
The cavity is assumed to be empty (no internal sources) and is illuminated 
by a plane wave in the 4> = 0 plane. For implementing the aforementioned 
solution, the cavity is subdivided into rectangular bricks and the results of the 
solution are shown in Figure 3. These are radar cross section (RCS) patterns 
and refer to a 1A deep cavity whose aperture is also 1A x 1A. The RCS 
pattern in Figure 3(b) applies to the cavity which is loaded with a resistive 
sheet of 100fl/D placed at its aperture, whereas the result in Figure 3(a) is for 
the untreated empty cavity. The simulation of the resistive sheet was done 
through direct discretization of the first order variation of the functional 
F as given in (7) and alternatively by modeling the resistive sheet as a 
dielectric layer of thickness A/20 having the dielectric constant computed 
from (2). As shown in Figure 3(b) the results based on the two simulations 
are in reasonable agreement and the differences among them is due to the 
finite thickness which was necessarily introduced in the physical model of the 
resistive sheet. As noted in [9] and [10], the dielectric layer introduces vertical 
components of the electric field which are not present in the resistive sheet. 
It is certainly of interest to point out that the presence of the resistive sheet 
at the aperture surface reduced the RCS of the cavity by lOdB at normal 
incidence and by as much as 20dB at grazing incidence. 

The second geometry which was considered is a circular metal-backed 
cavity again situated in a ground plane. In this case the aperture of the 
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circular cavity is loaded with a sheet having non-uniform resistivity given by 


R(p) = 


o.i + 10 (^) 


t oo 


p > 0.5a 
elsewhere on S a 


( 9 ) 


where a denotes the radius of the aperture. Results with and without resistive 
loading for a cavity having a = lin. and a depth of 0.25in. are shown in 
Figure 4. These RCS patterns were computed at 16GHz and the incident 
field was a plane wave polarized along the 0 or <f> directions. Again, the data 
in Figure 4 demonstrate the validity of the proposed mathematical model. 
Also, as in the case of the rectangular cavity the presence of the resistive 
cards reduces the overall RCS of the cavity and this reduction is primarily 
due to the rer^ed field intensity near its perimeter. 


4 Conclusions 

A formulation was derived for modeling resistive cards within the context of 
the finite element method. Essentially, the pertinent variational functional 
was supplemented with an additional boundary integral over the surface of 
the resistive sheet/card. Results based on the discretization of the functional 
were also presented and these were aimed at demonstrating the accuracy of 
the derived mathematical model. 
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FIGURE CAPTIONS 


Fig. 1 Cross Section of a dielectric volume enclosing a resistive card (a) card within the 
dielectric (b) card on the surface of the dielectric 


Fig. 2 Geometry of a cavity-backed aperture in a ground plane. 


Fig. 3 Monostatic RCS for a rectangular aperture (IX x IX) backed by a rectangular cavity 

(IX x IX x IX) in the <J>=0 plane, (a) Empty cavity without resistive loading (b) 
Empty cavity whose aperture is covered with a resistive sheet having a resistivity of 

lOOQ/ta. (Solid and dashed lines correspond to results based on the mathematical 
simulation; circles and squares refer to results based on the physical modeling of 
the resistive sheet) 


Fig. 4 Monostatic RCS for a circular aperture (1 inch in diameter) backed by a circular 

cavity (1 inch in diameter and 0.25 inches deep) at 16 GHz. (a) Empty cavity, (b) 
Empty cavity covered with a resistive sheet whose resistivity is given by equ. (a). 
(Solid and dashed lines correspond to results based on the mathematical 
simulation; circles and squares refer to results based on the physical modeling of 
the resistive sheet) 



Figure 2 

















